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Abstract. - We present an analysis of a system of weakly coupled Hubbard chains based on 
combining an exact study of spectral functions of the uncoupled chain system with a renormal- 
ization group method for the coupled chains. For low values of the onsite repulsion U and of the 
doping 5, the leading instability is towards a superconducting state. The process includes ex- 
cited states above a small correlation pseudogap. Similar features appear in extended Hubbard 
models in the vicinity of commensurate fillings. Our theoretical predictions are consistent with 
the phase diagram observed in the (TMTTF^X and (TMTSF^X series of organic compounds. 



Different classes of materials such as (i) high-T c superconductors, (ii) quasi-one dimensional 
(ID) organic conductors, and (iii) heavy-Fermion compounds show superconductivity (SC) 
phases whose physics is far from being well understood [1] . All these materials have as common 
feature the existence of non-Fermi liquid (NFL) metallic phases. In addition, in classes (i) 
and (ii) electronic transport is quite anisotropic, space dimensionality being expected to play 
an important role in the observed exotic spectral properties. On the other hand, SC phases 
were not observed in incommensurate quasi- ID materials [2,3]. 

The Hubbard model has attracted much interest as the simplest model for the descrip- 
tion of the electronic properties of these low- dimensional materials. While the ID Hubbard 
model has no SC phases for on-site repulsion U > [4], there has been evidence that the 
two-dimensional (2D) Hubbard model has a d-wave pairing instability [5]. In addition, the 
undoped ID Hubbard model is the prototype used in the study of the confinement problem [6] . 
Therefore, the study of the interplay of the crossover from one- to higher-dimensional systems 
with the electronic correlations is of great importance for the understanding of the exotic 
physics detected in the above materials [1]. An interesting toy model for that important prob- 
lem is a system of N c h — > oo weakly coupled Hubbard chains, each chain having N a —> oo sites 
and the hopping integrals such that fn > t±. Here ty and t± correspond to electron transfer 
along each chain and between first-neighbor chains, respectively. Importantly, we find below 
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that for electronic doping concentrations 8 ^ 0, where 5 = (1 — n) and n is the electronic 
density, the singlet SC spectral functions of the one-chain problem show singularities for en- 
ergy values just above a correlation pseudogap, E u hb, which equals the Mott Hubbard gap 
Emh for 6 — * [7]. Note that half filling for the coupled-chain system could correspond to 
a small effective doping S for the single chains, due to charge fluctuations introduced by the 
inter-chain hopping [8]. Indeed, the key point is that for low values of U/t\\ and of doping 5 
there is a regime where E u hb ~ \$\ 27rt|| + [8^/£/i||/7r] e^ 27Tt ^ u is small compared to t|| and one 
can choose tj_ such that E u hb < tj_ < tu. In such a regime the above finite-energy singularities 
of the one-chain problem give rise to relevant instabilities in the system of weakly coupled 
chains. Therefore, a complete study of the instabilities of the coupled chains having as start- 
ing point the properties of the one-chain problem must take into account both the low- and 
finite-energy singularities in the one- and two-electron spectral functions of the ID problem. 

In this Letter we combine the pseudofermion dynamical theory (PDT) for these spectral 
functions [3,9,10] with a renormalization group (RG) analysis [4,11,12] to study the insta- 
bilities of the above system of weakly coupled Hubbard chains. For the (TMTTF^X and 
(TMTSF)2X series of organic compounds, whose phase diagram is represented in Fig. 1, 
we identify E u hb with the experimental correlation gap [13] of Fig. 3 of Ref. [8]. Above 
the energies/temperatures of the inter-chain hopping and electron-phonon interactions, the 
ID electronic correlations are expected to dominate. Indeed, the ID Hubbard model has a 
NFL metallic phase for 5 ^ and a Mott-Hubbard insulator phase for S — 0, consistently 
with the occurrence of NFL and Mott-CO phases in the diagram of Fig. 1. This model was 
used successfully in the study of the photoemission spectrum of the incommensurate organic 
compound TTF-TCNQ [2,3, 10, 14]. The description of the above quarter-filled compounds 
involves extended Hubbard chains. However, due to dimerization [13] the low-T problem can 
be mapped onto our toy model for half filling, with U being an effective interaction parameter 
whose value depends on the bare on-site and first-neighbor (V) repulsions. Importantly, for 
energy and/or T values just above the dimerization gap, this leads to a correct description 
also for the small-<5 metallic phase, provided that in our spectral-function expressions 6 is 
replaced by an effective doping 5 = (1/2 — n) relative to quarter filling. (For instance, for en- 
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ergy /frequency (uj — Emh) > small, the optical conductivity reads Re o(uj) oc (uj — Emh) 1 ^ 2 
for both the 5 = ID Hubbard model [15] and the 5 = spinless fermion model with first- 
neighbor repulsion V [16].) We include in the perturbative expansion the inter-chain hopping 
t±. In this respect our method differs from alternative approaches where the expansion pa- 
rameter used is the electron-electron interaction U [17]. The latter perturbation scheme on 
U misses the effects that we discuss here. The influence of the inter-chain hopping terms has 
been studied by an expansion in the inverse of the chain coordination number [18]. 

The single Hubbard chain in a chemical potential \i reads H = H,so(4) — ^[N a — N], where 
^50(4) = t + UY.ftiA - V2] [n jd - 1/2], t = -t,| £ ji<T [c} CT c j+1(7 + h.c], N = J2 a N a , 
N a = J2j n j,<?> c jcr an d c ja are electron operators of spin projection a at site j = l,...,N a , 
and rij ia = (cj CT Cj CT ). #so(4) commutes with the six generators of the ry-spin and spin algebras 
[7,9,19] and hence the 77-spin value n and rj z = —8N a /2 are good quantum numbers. We 
consider the < S < 1 metallic and 5 = insulator phases with zero magnetization for the 
initial ground state (GS). The electronic charge reads — e and 2fcp = (1 — S)ir. We study the 
spectral functions B ( j\k,uj) where I = ±1 and B ( ^ 1 \k,uj) = \{f\d\{k)\GS)\ 2 5{u - uj ffl ) 

and B^~ 1 \k,uj) = J2f> \{f'\0-d(k)\GS)\ 2 5(uj +w//,o) for uj > and for uj < 0, respectively, 
whose singularities we detect and characterize. The / and /' summations run over the possible 
excited states and w/ ; o and ujf.o ar e the energies relative to the initial GS. We find that the 
more divergent singular features occur for onc-clcctron # = Ip, transverse spin density wave 
(SDW) ?9 = sdw, and on-site and extended singlet SC # = oss and i? = ess, respectively. 
The corresponding operators read 6 lp (k) = c fc;(T , O sdw (k) = J^k' c l+k'.i c k',h ®oss(k) = 

J2k> c fe-fe',T c fe',i' and Oess(k) = J2k> cos(k')[ck-k' ,iCk> ,1 - Ck-k',iCk',l]- Here c\ a and c fei<7 are 
spin d =t, I electron operators of momentum k. For the above ID model, electrons are related 
to rotated electrons by a unitary transformation and rotated- electron double occupancy D r is 
a good quantum number for all values of the on-site repulsion U [7,9]. For < <5 < 1, D r = 
for the GS and the first and second upper Hubbard bands (UHBs) are spanned by D r = 1 
and 2 excited states, respectively [7,9,10], whose lower limits are at energy D r E u hb- The 
correlation pseudogap E u hb is called E u in Ref. [10] and is such that E u hb = 4t|| sin(7r<5/2) for 
U/t\\ ->■ 0, E uhb = U - 4t|| cos(tt(5) for U » t\\, E uhb = U + it\\ for 5^1, and S uftb = E M h 
for (5^0. 

Our study of the above spectral functions B$\k,uj) takes into account all microscopic 
processes described by the exact PDT introduced in Ref. [9]. The excited states can be 
described in terms of occupancy configurations of holons, spinons, and c pscudofcrmions [9,10]. 
We denote the r\ spin (and spin) and ?7-spin (and spin) projection of these objects by s c and a c 
(and s s and a s ), respectively. The holons (and spinons) have s c = 1/2, a c = ±1/2, charge ±2e, 
and s s = (s s = 1/2, a s = ±1/2, and no charge degrees of freedom). The c pseudofcrmions 
carry charge — e and have no spin and 7/-spin degrees of freedom. We use the notation ±1/2 
holons (and ±1/2 spinons) according to the values of a c (and <r s ). The rotated-electron 
double occupation D r equals the number of —1/2 holons. In the vicinity of the singular 
spectral features the spectral functions can be expressed as a power-law expansion whose small 
parameter for a given value of the momentum is the energy deviation from the singular feature. 
For the one- and two-electron spectral functions only specific types of AD r = and/or AD r = 
1 processes contribute to the leading-order term of such a power-law expansion. (Here AD r = 
D r is the deviation from the ground-state value D r = 0.) We call these processes dominant 
processes. All remaining processes give rise to higher-order power-law contributions, controlled 
by exponents larger than those studied in this paper. The evaluation of the weight distribution 
in the vicinity of the singular features of the above two-electron spectral functions proceeds 
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as for the one-electron case reported in Ref. [10]. The dominant processes involve changes 
in the occupancies of the c pseudofcrmions, two-holon s c = composite cl pseudofcrmions, 
two-spinon s s = composite si pseudofermions, and localized and non-interacting —1/2 Yang 
holons, which are not part of composite pseudofermions [7,9]. There are no —1/2 Yang holons 
and cl pseudofcrmions in the D r = GS. The a — c, si, cl pseudofcrmions carry momentum 
q — q + Qa{q)/N a . Here q is the bare-momentum and Q a (q) is the phase-shift functional 
defined in Ref. [9], whose expression involves the bare- momentum two-pseudofermion phase 
shifts &aa'(q,q') in units of 2tt, where a, a' = c, sl,cl. (Such phase shifts are studied in 
Ref. [10].) The pseudofcrmions are related to the corresponding pseudoparticles of Ref. [7] 
by a unitary transformation [9] . The cl energy dispersion e cl (q) = E uhb + e° cl (q) such that 
\q\ < Stt plays an important role in our study and is plotted in Ref. [7]. 

Nearly the whole Cooper-pair addition spectral weight is generated by processes such that 
AD r — 0, 1, 2. These processes generate three separated parts B^^ (k,u>) for the SC spectral 
functions, B^~^ (k, uS) = ^2 D =0 x 2 B^^ {k,uj). At 5 = only D. r = 2 excited states are 

allowed and thus B < ^ 1 \k,u>) = £>^ + 2 ^ (fc, has no singular spectral features. For 5 > and 
U > the SC singular spectral features result from transitions to D r — 1 excited states. The 
most divergent singularities in the SC functions occur for d — oss, ess in B^^\k, to) and result 
from finite-energy transitions involving creation of a cl pscudofcrmion and corresponding to 
cj just above the two t' = ±1 singular branch lines E L ss (k) = e c i{q) where q — [l 1 2kp — k}. 
Such a process involves a r^-spin flip which transforms a +1/2 holon onto a —1/2 holon and 
leads to the deviations An = — 1 and Arj z = +1. By use of the method of Ref. [10], we 

find B [ ^ 1] {k,uj) = C { +l){k){uj - for low positive values of (w - E L s ' s (k)) where 

C^ + t V(fc) is real and positive. The exponent Css(k) = CiH^) = CssH - sucn that ( ss (k) = 
-l + E,=±iW K p/ 2 + ^c0ci(i2k F , q)} 2 +E,=±i{®sici(i>k F , q)} 2 refers to k e ([1-2%, tt) 
and involves the phase shifts defined in Ref. [10] and the parameter 1/2 < K p < 1 [20]. The 
multiplicative constant C^ + t V ' (k) vanishes for k = l'tt and in the limits 5 — > and U/t\\ — > 

for all values of k. While C { ^\l'-k) = for the oss and ess cl pseudofermion singular branch 
lines, creation of a —1/2 Yang holon leads in the metallic phase to a single S- function peak 
for — oss, B^ s !(t'7r,u>) — SN a S(u> — Euhb). This involves a rotation in ?]-spin space called 
■q-pairing mechanism [19] such that An — and An z = +1. This peak has a significant 
fraction of the oss first UHB weight but does not contribute to the coupled-chain instability. 
It is absent for the ess function, such that B^\ (l'tt, E^jj) = 0. Thus, for the latter function 
the whole first-UHB weight corresponds to the spectral feature generated by creation of the 
cl pseudofermion. Consistently although the oss and ess cl pseudofermion singularities are 
controlled by the same exponents, C^\,(k) > C^p ,(k). Therefore, in our analysis of the 
coupled chains we do not consider the oss singularity. 

The remaining more divergent one-chain singularities were already calculated in Ref. [4] 
and for small positive values of u read b[ ^ (t^f, — w) oc b[^ ^ (±kp, u) oc (oj)^ 1 " and 

Bi sdw(^ 2k Fi -w) oc B ( ^{±2k F ,u) oc (w) Csd ™. The exponents are given by ( sdw = -[1-K p ] 
and Ci P = -(3/2){l - (l/'3)[l/(2K p ) + K p /2}}. Here -1/2 < Csd w < and -1 < ( lp < -7/8, 
whereas ( ss (k) is a monotonous decreasing function of k such that ( ss (k) > for k = [1 — 2% 
and S > and ( ss (k) — > — 1 as k — > tt. For U > 0, Cip an d Css(k) are a decreasing and 
increasing function of S, respectively, and £ip = —7/8 and ( ss (k) = —1 for 6 —> 0. The regime 
of interest for our problem corresponds to small values of S and U/t\\ and of [tt — \k\] > 
for C ss (k) 7 when E L ss {k) « 2t|| [vr — |fc|] + E uhb . The correlation pseudogap E uhb is small for 
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Fig. 2 - The exponents £i p and Css([l — S/4]n) as a function of 5 for U = t\\ (full line) and U = 1.5 1\\ 
(dashed line). 



S e (0, 0.02) and U/t\\ £ (0, 2), when E uhb w 52-Kt\\ 



-27rt||/C7 <K t and the 



condition E u hb < ij_ < t\\ is fulfilled. The exponents ( ss (k) for k = (1 — S/4)tt and Cip are 
plotted in Fig. 2 as a function of <5 for U/t\\ = 1.0 and 1.5. (The exponent ( s du. > — 1/2 for 
{/ > does not appear in the figure.) For low values S below some value So = 5o(U/t\\) the 
exponent £ ss is the smallest and for S > So, Cip is the smallest exponent. Let k = (1 — S/a)ir 
with a > 1 in the argument of ( ss (k). For the same value of U/t\\ (and a), <5 is an increasing 
function of a (and U/t\\). Thus, the same qualitative results are reached for smaller values of 
U/t\\ and So, when the inequality E u hb « t± « t\\ is fulfilled. 

The renormalization of inter-chain couplings can be directly expressed in terms of the 
exponents which characterize the intra-chain spectral functions. For the ID system the most 
divergent singularities occur in the lp and ess functions and then we can restrict our attention 
to terms of the type, 



H jc = J /( fc/ )/( fc ") c I,fc-fe',T4fe'a c »',fc-fe",T c i',fc'M ; 



(1) 



Here Hj c and H p h refer to Josephson coupling and particle hopping, respectively, f(k') defines 
the extended s pairing, and i and i' denote chain indices. These couplings are relevant and 
grow at a rate determined by the RG equations: dt±/dl = — d p t± and dJ/dl = —(ss J+Ct\. 
Here I = — log(A/A ), A is the initial value of the cutoff, t± and J are the couplings divided 
by A, and C is a constant of order unity. Initially, the rate of growth of the Josephson coupling 
is proportional to t\, and when this coupling reaches a value, in dimensionless units, of order 
its growth will be determined by the scaling dimensions. The analysis is valid until one 
of the couplings, tj_ or J, becomes of order unity, or when A ~ E u hb- Integrating the RG 
flow, we obtain i±(l) = t ± (0)e^ 1 and J (I) = P ± (0) [e 2 ^ 1 - e^ l }/[2\( lp \ - \( ss \] + J e^ 1 , 
where J is the initial Josephson coupling. Note that as for spin singlet SC, all four-operator 
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Hamiltonian terms have couplings whose leading-order term is tj_(0) [e 2 ' ( = lr ''']/[2|Cip| — 
However, since is largest for d — ss, the multiplicative constant £j_(0)/[2|Ci p | — |Ctf|] 
is also largest for i9 = ss. This is consistent with our choice of Hamiltonian terms, given 
in Eq. (JJJ. As 2|£i p > £ ss |, the RG analysis breaks down at a scale log(A e ///Ao) ~ 
— |Cip| _1 log[fj_(0)]. Provided that A e // is greater than E u m>, the analysis presented here 
gives a reasonable description of the properties of the system at energies higher than A e / / . 
Otherwise, the scaling equations cease to be valid at A ~ E u hb- Let us assume that A e // 3> 
Euht,- At this scale both t± and J become of order unity. Thus, the system presents two 
regimes: (i) at high energies or temperatures, greater than A e //, the system is best described 
in terms of decoupled chains: the effective inter-chain coupling, in absolute units, decreases 
upon lowering the temperature or energy, as ij_(T) ~ ij_(Ao) (T/Ao) 1- K Ip ', intcr-chain pair 
fluctuations increase, and the coherent contribution to the conductivity decreases; (ii) At 
sufficiently low energies or temperatures, as compared to A e //, the leading terms in the 
effective Hamiltonian are a one-electron and a pair-hopping term of similar magnitude, ~ A e jy. 

Note that the fact that the Josephson coupling has the largest scaling exponent ensures that 
its value at the point where the RG breaks down is the largest among all possible interchain 
interactions. The ground state is a strongly coupled superconductor which will show deviations 
from the BCS theory. The low-T phase can be viewed as realization of the inter-layer pair 
hopping model derived from coupled systems with strong RVB correlations [21]. Our analysis 
helps to delimit the region of validity of this hypothesis: (I) t± <C tn; (II) S or 5 close to 
zero ; (III) E u hb should be smaller than t±. Note that the states which contribute to the 
pair- hopping instability change the r\ value, i.e. the initial GS and the excited states are not 
related by a mere rotation in 7/-spin space. Thus, the present pairing mechanism cannot be 
considered directly related to the processes analyzed in the SO (5) theory of the formation of 
d-wave SC in the cuprates [22]. The study of Ref. [12] also leads to singlet SC, but does not use 
the singularities of the uncoupled system. Our results are universal in the sense that from the 
RG scheme used here we can infer that the low-T phase of the coupled-chain system will show 
long-range SC order. However, the precise nature of this phase, and the symmetry of the order 
parameter is dependent on the arrangement of the chains within the material. Our results 
confirm that in ID NFLs weak inter-chain hopping can induce SC and are consistent with the 
phase diagram of the series of quasi-lD compounds of Fig. 1. At low T these materials show 
a spin Peierls (SP) or SDW phase. Under pressure the (TMTSF^X compounds are driven to 
a SC phase, which is removed again if one further increases the pressure. Consistently, for low 
values of 5 we find that E u hb > t± away from the small U/tu region and according to the RG 
scheme of Ref. [4] the leading instability is towards a SDW state. Increasing pressure means 
here increasing tu and t± without changing both U and the condition t±/tu << 1. This drives 
the system into a small-C/ /tn region, where E u hb becomes smaller and t± is larger, and thus 

Euhb < t±, and the leading instability is towards a SC state. Since in the above B^}\k,uj) 
expression cits\k) decreases with decreasing U/tu and vanishes as U/tu — > 0, further increase 
of pressure removes this state. In the absence of pair hopping, the renormalization of the inter- 
chain conductance ceases at A e //, leading to an effective Fermi- liquid (FL) description at low 
T. The ability of our toy model to describe the diagram of Fig. 1 for X — CIO4 can be further- 
checked by deriving the values for its effective parameters from the NFL phase for T w 10 K 
and verifying whether they are consistent with the occurrence of the SC phase. The use of 
the PDT for the uncoupled system leads to a w-dependent conductivity for the NFL phase 
including a small lj — Drude peak and a finite-energy absorption for E opt < uj < E opt + 8tu 
where E opt = E u hb + e?i(0). This absorption is generated by creation of a —1/2 holon and 
two c pseudofermion holes at q and — q for |g| £ (0, 2kp w 7r). As found in Ref. [15] by use of 
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a preliminary version of the PDT, quantitative agreement with the X — CIO4 and T = 10 K 
curve of Fig. 1(C) of Ref. [8] is obtained provided that t\\ ~ 0.125 eV, U/tn w 1.5, and 
5 « 0.005. Importantly, according to Fig. 2, these values are consistent with the occurrence of 
the SC phase for smaller T. Tight-binding model (TBM) calculations lead to t\\ w 0.250 eV [8] 
because the width of the c energy-band filled sea for \q\ < 2fep ~ tt is about twice that of 
the naive TBM band [7]. Our mechanism also explains the absence of SC phases in quasi- ID 
materials whose 5 value is not close to commensurability. For TTF-TCNQ, 6 = 0.41 [2,3,10], 
the correlation pseudogap E u hb is finite for all values of U , and thus the condition E u hb < t± 
cannot be fulfilled. Although our analysis refers to tu > t±, our results suggest that upon 
increasing ixA|| the SC instability will be extended to larger values of both 5 and U. While 
for the quasi-lD system the instability corresponds to momentum values in the vicinity of it, 
we expect that for the high-T c superconductors square lattice such an instability will appear 
in the vicinity of [it, tt] as t±/tu — > 1. 

We thank C. Bourbonnais (C.B.), A. H. Castro Neto, R. Claessen, and L. Degiorgi for 
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